Strong-coupling scenario of a metamagnetic transition 
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We investigate the periodic Anderson model in the presence of an external magnetic field, using 
dynamical mean-field theory in combination with the modified perturbation theory. A metamagnetic 
transition is observed which exhibits a massive change in the electronic properties. These are 
discussed in terms of the quasiparticle weight and densities of states. The results are compared with 
the experimental results of the metamagnetic transition in CeRu2Si2. 



By "metamagnetic transition" , we describe an anoma- 
lous behavior of the magnetization as function of external 
field 6ext, namely a sudden increase at a finite field 6*^^. 
Many rare-earth materials exhibit this kind of behavior. 
One has to distinguish between those materials which al- 
ready show long-range (e.g. antiferromagnetic) order for 
zero-field and those which are paramagnetic. 

In the first case, the materials already have finite lo- 
cal moments at zero-field. Here a metamagnetic transi- 
tion occurs when the external field is stronger than the 
internal (antiferromagnetic) exchange between these mo- 
ments. This situation is found for example in CeFe2- 
based alloy sld. 

The other possibility, where a paramagnetic substance 
enters a high-magnetization state at a critical field &*xt 
is realizecLbiy.|aome heavjfcjfermion metals as for example 
CeRusSiaBM or UCoAffl. As will be discussed further 
below, experiments indicate a substantial change in the 
electronic structure at the transition, which is not yet 
fully understood. 

In this paper, we investigate a similar transition found 
in a relatively simple electronic model of heavy-fermion 
compounds. We examine the periodic Anderson model 
(PAM) with an external magnetic field. The periodic 
Anderson model describes the interplay between strongly 
correlated localized electrons with a band of uncorrelated 
conduction electrons. These two electronic sub-systems 
are coupled by a hybridization term. The Hamiltonian 
of the PAM reads: 
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on-site Coulomb interaction strength between two /- 
electrons. Throughout this paper, the conduction band 
will be described by a free (Bloch) density of states, 

Po{E) = jiJ^k^i^ ~ ^(^)): of semi-elliptic shape. Its 
width W — 1 sets the energy scale, and its center of 

gravity the energy-zero: Tu — ^(^) ^ 0. The mag- 

netic field bcxt is given in energy units of the band width 
W i z„ = +!(-!) for <7 ==T (i))- The external field 
couples equally to the spin of the / and conduction band 
electrons. This is in our opinion most appropriate for the 
model Hamiltonian (|l|) where the /-states are taken to be 
non-degenerate (s-type). Without using the full orbital 
degeneracy, one could alternatively use the g-factors of 
the real materials. Other authors have even completely 
neglected the coupling of the magnetic field to conduc- 
tion band states arguing that the corresponding g factor 
is negligiblelZloEl. 

We e*aploy the dynamical mean-field theory 
(DMFT)ll£l in combination with the modified per- 
turbation theory (MPT)E2l to determine the one-electron 
Green's function, from which the excitation spectrum as 
well as magnetization, effective mass and other quan- 
tities can be calculated. This inethod has previously 
been applied, to the paramagneticll3 and the ferromag- 
netic PAMI13, so we can confine ourselves to a short 
summary: The underlying idea of the DMFT is that the 
local self-energy suctuas occurs in the limit of infinite 
spatial dimensionaljll3, can be taken to be that of an 
appropriately defined single-impurity Anderson model. 
^-We solve the latter using the modified perturbation 
theory. Orp-*tarting point is the following ansatz for the 
self-energyllaO: 
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Here, sg^ (fia) and st (fj^) are the creation and anni- 
hilation operators for a conduction electron with Bloch 
vector k and spin a (a localized electron on site i and 
spin a) and n'^ = flf., (sj:^ = ^ Efc e^'^^'s..). The 

dispersion of the conduction band is e(A:) and ef is the po- 
sition of the localized level. The hybridization strength 
V is taken to be fc-independent, and finally U is the 



Oct and Pa- are introduced as parameters to be determined 
later. Ect^*^*^'' (E) is the second-order contribution j-to 
perturbation theory around the Hartree-Fock solutionta. 
Equation (^ can be understood as the simplest possible 
ansatz which can, on the one hand, reproduce the per- 
turbational result in the limit J7 — *■ 0, and, on the other 
hand, recovers the atomic limit for appropriately chosen 
a a- and /jjla. 
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Using the perturbation theory around the Hartree- 
Fock solution introduces an ambiguity into the calcula- 
tion. Within the self-consistent Hartree-Fock calculation, 
one can either choose the chemical potential to be equiva- 
lent to the chemical potential of the full MPT calculation, 
or take it as parameter fl to be fitted to another physi- 
cally motivated constraint. In reference |l^ the Li^ mger 
theoremtj, or equivalently the Friedel sum rulcEHlEiJ, was 
used to determine fl. As discussed in Ref. ^ we use 
the physically motivated condition of identical impurity 
occupation numbers for the Hartree-Fock and the full 
calculation {n'^J'^^^ — n'J'^) to determine /i, which also 
allows for a consistent extension of the method to finite 
temperaturealS'Ed. Except for symmetric parameters this 
will lead to ap,approximate fulfillment of the Luttinger 
theorem onlyllil. The key features of the results pre- 
sented below, however, are not decisively influenced by 
this shortcoming. This was checked by adopting the al- 
ternative approach of fitting fl using the Friedel sum rule. 
A more detailed analysis of the different possibilities to 
determine fl is found in reference 22 where the DMFT- 
MPT was applied to the single-band Hubbard model. Fi- 
nally, the parameters and 13^ have to be determined. 
Instead of using the "atomic" limit of ^ = as was done 
for example in references 16 1^,^, we make use of the 
moments of the spectral density. This procedure is de- 
scribed in detail in references |ll|,^ The result not only 
fulfills the = limit, but also recovers the high-energy 
behavior of the Green's function up to the order (■^). 

It has been shown that the approximation scheme, as 
described above, gives qualitatively reliable results by 
comparing with numerical renormalization group the- 
ory and quantum Monte Carlo calculations for two 
different "strong-cauoling effects", namely the Mott- 
Hubbard insulatorEHlcipid ferromagnetism in the single- 
band Hubbard modcSEa. The periodic Anderson model 
was investigated-, by this and related methods-ip, the 
paramagnetidlaEZi and the ferromagnetic phasetaEa. In 
the following we present results for the periodic Ander- 
son model in the paramagnetic, but ^close to the ferro- 
magnetic region of the phase diagranJlj. 

In Fig. |l], the f-magnetization m^^'' — - 



(t) 



(f) 



and the 
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(spin-dependent) quasiparticle weight Zc 

— ^"^g^^ l-E=o)~"^ are plotted as function of the exter- 
nal field &ext for ef = -0.4, [/ = 4, n(*°*) = 1.85, V = 0.2 
and zero temperature. These parameters away from the 
symmetric point (where ef = — -j) are close to the ferro- 
magnetic phase whichrwould be reached for lower electron 
density (n(*°') < 1.7)ll^. The increase of m(&oxt) (thick 
solid line) shows a sharp rise at &oxt = &oxt ~ 0.033. 
We call this phenomenon a metamagnetic transition. In 
Fig. 1^, a discontinuous jump in to^^' is visible. However, 
due to the finite numerical resolution of our calculations 
we cannot rule out a continuous transition. For a range of 
^oxt < ^oxt a high-magnetization and a low-magnetization 
self-consistent solution co-exist. For T = 0, the latter is 
always stable. The spin-| quasiparticle weight vanishes 
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FIG. 1: m'''(6cxt) (solid line) and the (spin-dependent) quasi- 
particle weight Za as function of &ext for ef = —0.4, [7 = 4, 
^(tot) ^ ^ y ^ Q 2 and r = 0. The inset shows n«^' as 
function of 6cxt. 
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FIG. 2: Same as Fig. |l|for 6cxt = 0.1 as function of U. 



at implying a maximum in the effective mass. For 
&oxt > ^oxti the quasiparticle weight is significantly larger 
than for &oxt = 0. As shown in the inset, the metamag- 
netic transition is also accompanied by a sharp increase 
in /-occupancy n'^^ at foext = ^ext- 

In Fig. H, the magnetization and the quasiparticle 
weights are shown as function of interaction strength U 
for fixed 6oxt = 0.1. At U — U* a metamagnetic transi- 
tion is also observable, being accompanied by the same 
behavior of Za- as discussed above. 

The /-DOS for 6cxt £ {0.0,0.02,0.04,0.06} is plot- 
ted in Fig. H for the same parameters used in Fig. [|. 
For zero-field, the DOS consists of lower charge excita- 
tion at ef, upper charge excitation at ef -I- ?7 and the 
Kondo resopance at E ^ fi which is split by the co- 
herence gapEHl. Applying a small external magnetic field 
(bcKt < 0.04) induces a Zeeman shift proportional to 6oxt, 
which is best visible in the upper charge excitation. The 
Kondo resonances in the spin-f and | channel seem also 
to be shifted, however, in the opposite direction: Spin- 
t to higher, and spin-| DOS to lower energies ("inverse 
Zeeman shift"). This is very unexpected behavior and 
is only found for parameters that show a metamagnetic 
transition. For other parameters further away from the 
ferromagnetic phase (e.g. U < U* or also at the symmet- 
ric point Cf = — a "normal" Zeeman shift is observed 
for the Kondo resonancecj. Closer investigation shows 
that the apparent shift is indeed no shift, but a suppres- 
sion of spectral weight in the spin-f channel just below 
the Fermi energy. This suppression is also indicated by 
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FIG. 3: /-DOS for the parameters of Fig. |l| with boxt as in- 
dicated in the different rows. The solid (dotted) line corre- 
sponds to the spin-| (l) channel. The two columns show the 
same DOS at different energy ranges: lower charge excitation 
and Kondo resonance in the left and the upper charge excita- 
tion in the right column. The chemical potential is given by 
the arrows. 

a decrease of Z-^ as seen in Fig. ^ 

For fields above the metamagnetic transition, the pic- 
ture changes dramatically. The lower charge excitation is 
fully saturated and the shift between the two spin chan- 
nels of the upper charge excitation is much larger than 
one would expect from an extrapolation from the other 
figures. Finally, the Kondo resonance disappears in the 
spin-f DOS. A broad structure located well above the 
chemical potential is visible in the |-channel only. Thesa 
DOS strongly resemble those of the ferromagnetic PAMO 
whereas the DOS for lower-ields more or less correspond 
to the paramagnetic DOSO plus the Zeeman shift. So 
the metamagnetic transition, closely resembles the para- 
to ferromagnetic transitionll^. This agrees with the fact 
that the metamagnetic transition is clearly a strong cou- 
pling phenomenon as was also found for ferromagnetism 
in the intermediate- valence regimecS (see Fig. 

Let us shortly discuss the nature of the metamag- 
netic transition exhibited by our results. The origin 
of the transition is somewhat different from previously 
discussed approaches to metaniaaaf^m, such as the 
Kondo-volume collapse approachEJly'L^I, or exnlajTtaiioiis 
based on special features of the density of stateaafla'tSE^. 
For example, solviag-the PAM within self-consistent per- 
turbation theoryESEa, one also finds a metamagnetic 
transition. Here, however, the transition originates from 
sharp features in the excitation spectrum which are al- 
ready present for zero-field. On applying the field, the 
Zeeman shift can push these features across the chem- 
ical potential and thus dramatically change the spin- 
dependent f-occupation and hence the magnetization. In 
our results however, the transition can be traced back 
to dramatic changes of the correlation-induced features 



of the excitation spectrum due to the external magnetic 
field. This clearly distinguishes the scenario presented 
in this paper from previously discussed metamagnetic 
transitions in the PAM. This manifests itself furthermore 
in the observed hysteresis which cannot occur in a sce- 
nario where the transition is due to the Zeeman shift of 
a strongly peaked DOS. 

A phcnomenologically similar metamagnetic trau.sit.ioji 
is known to exist for the half-filled Hubbard modelE3'LElEZl. 
Here, a jump in the magnetization as function of exter- 
nal field is found for interaction strengths U close to 
the critical Uc separating the Mott-insulator from the 
metallic regime. Similar to the transition described in 
this paper, the transition in the Hubbard model also 
shows a hysteresis and is indicated by a suppression of 
the Kondo resonance. Another similarity might be seen 
in the fact that in both metamagnetic transitions, the 
high-magnetization state is characterized by a stronger 
localization of the correlated electrons. In the Hubbard 
model, this transition is found for half-filling where it 
leads from a metallic to an insulating state. In the tran- 
sition presented here we also see a tendency towards 
stronger localization of the f-electrons. Contrary to the 
situation in the half-filled Hubbard model, however, the 
PAM is metallic above and below the critical field. Fur- 
thermore, the suppression of the Kondo resonance occurs 
only in the proximity of the ferromagnetic phase, and 
only in the spin-l channel. 

The metamagnetic transition is not simply an enlarge- 
ment of the 'ferromagnetic' phase caused by the mag- 
netic field. If this were the case one would expect simi- 
lar behavior using other approximation methods, which 
yield |almost the same ferromagnetic phase diagram as the 
MPTca such as the modified alloy analpgsi-aBjd the spec- 
tral density approximation for the PAMMESo. However, 
both methods do not show the metamagnetic transition. 
The low-energy ("Kondo" Ix)hysics, which are not recov- 
ered by the other method£j play a decisive role. 

What is the relevance of our results to the metam- 
agnetic transitions observed experimentally? The most 
well-known example of a heavy-fermion compound show- 
ing a metamagnetic transitiejaJrom para- to a "ferromag- 
netic" state is CeRu2Si2EJE3cll. This material exhibits 
peculiar behavior alongside the metamagnetic transition. 
The effective mass shows a sharp maximimi at 6oxt = Kxt^ 
and for 6oxt > ^cxt suppressedaaati, which is ex- 
actly the behavior of the metamagnetic transition of the 
PAM discussed above. Furthermore, the experimentally 
conffrmcdcl stronger localization of the /-electrons in the 
high-magnetization state is indicated by the jump in n^^^ 
that we discussed above.. 

There are two findings for CeRu2Si2, which are not 
reproduced pby our approximation. Neutron scattering 
experimentaij revealed antiferromagnetic inter-site cor- 
relations for &ext < ^ext- These correlations can not be 
found within our approximation scheme since it is based 
on dynamical mean-field theory. However, the param- 
eters we have used above are close to the antiferromag- 
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netic regime ofthePAMea. It seems therefore reasonable 
to expect antiferromagnetic correlations in the zero-field 
(low-magnetization) state. The second ^flfect is the vol- 
ume effect of the metamagnetic transitionOcS. A Kondo- 
volume collapse similar to that discussed in connection 
with the 7 ^ a transition of CeEj has been put for- 
ward as possible source of the metamagnetic transition 
in CeRuoSio PHp^l Our model does not include any cou- 
pling to the lattice. However, the significant change in /- 
occupancy that we find should result in a volume change, 
if the lattice coupling were to be included. This could 
even amplify the transition. 

Contrary to our results, the metamstgnetic transition 
in CeRu2Si2 does not show hysteresisQ. Let us point 
out that in our results wc find merely a "mathemati- 
cal hysteresis" . In Fig |l|, the thick line is the physically 
meaningful result whereas the thin line represents only 
a mathematical solution to the equations to which we 
cannot ascribe a clear physical meaning. However, one 
could assume that for finite temperatures, a true hys- 
teresis would be found. Experimentally, a metamagnetic 
transition from a paramagnetic state ttat is accompa- 
nied by a hysteresis is realized in UCoAIeI. This material, 
however, does not show a shapp maximum in the effec- 
tive mass of the quasiparticleaa so it is unclear whether 
the picture described in this paper is relevant for this 
material. 

At this point, let us comment on the critical fields nec- 
essary to drive the metamagnetic transition. In our units, 
bext = 0.033 corresponds to several hundreds of Tesla if 



the bandwidth would be of order leV. This is much 
too large compared to the above-cited experimental re- 
sults. We believe this is a limitation of the MPT ap- 
proximation. In our explanation, the transition is closely 
connected with the low-energy properties, the existence 
and the width of the Kondo resonance. The MPT-iiow- 
ever, tends to overestimate the low-energy scaleai3, so 
the absolute values of the critical field as obtained by 
our calculation should not be taken too seriously. 

To summarize, we have presented a new scenario for 
a metamagnetic transition based on a strong electron- 
electron coupling effect. For a periodic Anderson model, 
with parameters located in the paramagnetic regime be- 
tween the ferro- and the antiferromagnetic phase, an in- 
crease of the external field for sufficiently high interaction 
strength leads to a sudden sharp increase of the magne- 
tization (metamagnetic transition). This is accompanied 
by a suppression of low-energy spectral weight (Kondo 
resonance) in the spin-t DOS. The features of the tran- 
sition show strong similarities to the experimentally ob- 
served metamagnetic transition of CeRu2Si2. The de- 
creasing spectral weight near the Fermi level prior to the 
metamagnetic transition could serve as an experimen- 
tally accessible indicator for the relevance of the proposed 
mechanism to the real physics of CeRu2Si2. 
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